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Abstract 


The effective level of line formation for spectroscopic absorption lines 
has long been regarded as a useful parameter for determining average atmo- 
spheric values of the quantities Involved in line formation. The identity of 
this parameter has recently been disputed. Here we re-establish the depen- 
dence of this parameter on the average depth where photons are absorbed in a 
semi-infinite atmosphere and show that the mean depths derived by others are 
similar in nature and behavior. 
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The concept of a "level of line formation" is scarcely novel in 
Astrophysics (see the review by Stromgren, 1951, p. 204). Nevertheless, this 
concept was not introduced into studies of planetary physics until 15 years 
ago, in a theoretical paper concerned with the hazy atmosphere of Venus 
(Chamberlain, 1965, p. 1190). Recent ly, the usefulness and functional 
behavior of this parameter have been disputed (Wallace and Hunten, 1978; 
Kattawar, 1979). We contend that for a homogeneous, semi-infinite atmosphere 
the effective level for line formation is a straightforward matter flowing 
from the transfer equation. 

We define the mean depth as the first moment of the sink function 


Iq S(t) xdx 
Jo S(t) dx 


( 1 ) 


where 


S(x) = (1 - 5>) J(x, u Q ) , 


( 2 ) 


<*> is the single scattering albedo, and J is the mean intensity averaged over 
direction at depth x. In Chandrasekhar's (1950) first approximation for a 
semi-infinite atmosphere we have 

J(t, p 0 ) = (Le _kT + Ye'^o) (3) 


where L and y are functions of p^, p q , and £, ttF is the incident flux; 

= 1//T and k * (l - w)^^/pj. Substitution of (2) and (3) into (1) gives 
(adapted from Chamberlain, 1965) 
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<T> 


(4) 


' (1 -“i)I72 * "o' “l 1 

Note that the first term on the right hand side is the asymptotic limit 
for <t> for weak lines with a weakly absorbing continuum (i.e., « *1). 
However, without loss of generality, x can represent optical thickness in the 
continuum as well as a line. 

A physical interpretation of this result comes from considerations of a 

random walk of photons (Chamberlain, 1978, p. 140). Very crudely, the number 

2 

of scatterings required for a photon to migrate over depth x » 1 is x . With 
a probability of absorption at each encounter of (l - Si), we would expect to 
find a substantial fraction of the photons absorbed only when 

x 2 (l - a) • 1 

or x » (l - a)" 1 / 2 . 

Although the sink function S(x) = (l - a) J (x, y Q ) technically depends 
on y Q , in the limit of weak absorptions and non-grazing incidence, <x> is 
practically independent of y Q , because the y or "solar term" (representing 
primary scattering) in Equation (3) is relatively unimportant at the great 
depths of <x>. 

Wallace and Hunten (1978) derive the equivalent width W, which comes from 
wavelength integrals of intensity I v (x v = 0; y, y Q ), both in the line and 
continuum, and therefore depends on the viewing geometry. The problem lies in 
assigning a pressure or depth where the effective pressure-broadening 
parameter, a g , can be evaluated. In a real atmosphere a varies with depth, 
but we suppose there is an effective depth x g where the local «(x g ) = best 
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represents the entire composite line. The concept Is reasonable, but their 
evaluation of a Is suspect. They first of all equate equivalent widths for 
constant-pressure and varylng-pressure models In the limit of a vanishing 
albedo, 5 ♦ 0, or a single-scattering atmosphere. 

As noted above, a u, y Q dependence Is then built In; but the limit of 
single scattering scarcely yields a radiative-transfer (multiple scattering) 
model. For larger 5 they obtain the numerical coefficient of <* e by matching 
curves of growth for the two models (constant and varying pressure), but the 
directional geometry remains the same because it is tied to the case 5=0. 
Through a trial -and-error method, they multiply the two curves of growth by 
different factors which are functions of y, y Q , and 5 to achieve overlap and 
derive: 


j — _2_ f — . i- j/2 . i— i.—M 1 

* T C l (l - IFmT 1 ♦ U -!) J 


(5) 


where a/t c is the line width at unit optical depth. It is surprising to see a 
reflecting layer air mass factor n = y -1 + in a formula for a scattering 
layer. Chamberlain (1965) showed that for multiple scattering the y, y Q 
dependence for weak line equivalent widths, was (y + y Q ). However, if we 
assume that their method of matching curves of growth is correct, it is easy 
to suggest an alternative formula with the correct (y + y Q ) variation. The 
results and figures presented by Wallace and Hunten are for y = y Q ; in this 
case, the factors (y + y Q ) and l/(y~* + y Q ”^ ) are identical in variation . 
Therefore, since they use arbitrary numerical factors to force their constant 
and varying-pressure models to give identical curves of growth, we can 
multiply Equation (5) by (y -1 + + u Q )/4 which is unity when y = y Q . 


Their formula then reads 



ill - 
(1 - “) 


T7? 



( 6 ) 


Not only does this give exactly the same curve of growth as (5), when y * y Q , 
but it also contains the Chamberlain (1965) y ^-variation. 

Kattawar's (1979) quarrel with Chamberlain's <t> is conceptual. His own 
definition of the effective level of line formation is 


T eff C tN se < ,)dT ' C N se (T)dT 


(7) 


where N se (t) represents the number of photon scatterings occurring at depth t 
for light which emerges. In other words, is the average depth at which 
an emerging photon has scattered. On the other hand, Chamberlain's <t> 
represents the average depth at which photons are absorbed. 

Kattawar's principal argument against Chamberlain's <x> was that as <*> 
approaches unity, <t> goes to infinity. This is a result of the fact that 
as S> approaches unity the light intensity deep in a homogeneous medium becomes 
a constant with depth. Obviously, the depth-averaged value of such a constant 
function will tend to infinity. 

However, it is a simple matter to show that Kattawar's T e ff. as defined 
in Equation (7), has the same asymptotic behavior as <t>. Let N s (x) be the 
total number of scatterings occurring at depth t. N $e is less than N $ by a 
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factor P(£>, t) which represents the probability that a photon, which has 
scattered at depth t, will eventually random walk out of the atmosphere. 
Following the lead given by Chandrasekhar (1943) for the random walk of 
photons with an absorbing screen we obtain (see Appendix) 

P(&, x) = exp [- (t^ In (l/i£))l/2] (8) 

Multiplying this probability by the mean intensity, J(t,Pq), gives N ge 
(times a constant) for insertion into the integrands of Equation (7). We then 
have 


[k + m]~ 2 + ~ [u 0 _1 + m] -2 
[k + sn]~* + ^ [uq - 1 + m]" 1 


(9a) 


where m = [In (1/S i 2 )]* . Note that for u> * 1, ^ - ^)/v*g and 

1/2 

m ■* u^k(2) . Since k ♦ 0 as S> 1 , we have 


Teff = (i" + + ,J ° ‘ 


(9b) 


This approaches infinite depth as w 
Chamberlain's <t> in Equation (4). 

This raises the obvious question, 
Kattawar argued otherwise? The answer 
The basic ;r.T1_n is t h it his Mol 


ri«ars 1 only slightly slower than 


if I,, is so similar to <t>, why has 

CM 

is found in his calculat *' oal v J el. 
■ :i '0 is filii'e. ibis is ' ot 
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appropriate for a comparison with <t>, which Is derived for a semi -finite 
homogeneous medium. 

The Identity of is further obscured because the top tenuous cloud 
layers are weighted the same as the deep atmosphere since his Integration step 
size Is In kilometers (h e ff) and not optical depth. 

In Figure 1 we have compared a typical example of Kattawar's radiance 
distribution with a classic source function. In effect, these are the two 
weighting functions for the calculations of (A) and <t> (B) In Equations 
(7) and (1) respectively. The effects of the finite atmosphere on are 
evident. As » approaches unity, will remain finite. However, with the 
above defects corrected, Kattawar's would match <t> in behavior. 
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Appendix I 

The Escape Probability 


In order to calculate the probability of escape from a semi -Infinite 
cloud for a photon which at time t ■ 0 scatters at depth t, we refer to the 
classical random walk problem with an absorbing screen as put forth by 
Chandrasekhar (1943). Then, the probability of escape in precisely n 
scatterings is given by 

P («. t, n) • e‘ t2/2n . (10) 

Therefore, the total escape probability Is 


P(i, t) 


jo . 

/: * h &)‘ /2 «- ,J/2n «» 


(11) 


Since P(«* t, 0) ■ 0 (l.e. the probability for escape from depth t with zero 
scatterings Is negligible) then we can extend the integrals' lower limit to 
zero. Thus, we have for the numerator 

I„ • T r, n- 3/J e- A/n e-™ dn 

where A » t 2 /2 and m * -In (u>). Consider 

f(A) • /;e' A/n - nn n ' 3/2 dr. (I?) 


e yStfWfciV .. 
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Me see that 


f(A> • - If i; e’* /n - m *- m dn. 

Letting x ■ (n) 1 ^ 2 , we have 

f(A) . . 2 d. /* e-*/* 2 ’ ”“ 2 dx. 

Now letting x * (A)^ 2 y, we obtain 

f(A) • -2 ^ [(A) 1/2 g( tj ] 

where 

g(r) • i; e-v 2 - 1/yZ dy 

and Y * mA. Now, considering g{ir), let z « 1/y 2 . We then derive 

g( T) ■ i; ^ - * Z dz/z 2 

so that 

gtrl.-^e-^-^dz]. 

1 /2 

Finally making the change of variable z = (y) ' y, we have 



Substituting from Equation (14) then gives 

g(r) ■ - [(t) 1/2 g(t)) 

with a solution of 

g(v) » C-f 1/2 exp {-2y 1/2 > 

where C Is a constant of Integration. Substituting back Into Equations (13) 
and (12) gives 

P(®» t) - exp [- (t 2 In (1/G^)] 1 ^ 2 ]. 
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Figure Caption 


Figure 1. Optical depth dependence of (A) a radiance distribution for a 
finite atmosphere (Kattawar, 1979; figure 2b) and (B) a source function 
for a semi-infinite homogeneous atmosphere (Chamberlain, 1970). 
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